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THEORETICAL  SOLUTIONS  FOR  A  STEADY  LAMINAR  PLOW 
OP  A  COMPRESSIBLE  FLUID  IN  THE  ENTRANCE  REGION  OP  A  TUBE* 

by 

Tau-Yi  Toong**  and  Joseph  Kaye*** 


Massachusetts  Institute  of  Technology 


SUMMARY 

The  boundary-layer  equations  of  continuity,  momentum, 
and  energy  for  the  steady,  laminar  flow  of  a  compressible 
fluid  in  the  entrance  region  of  a  round  tube  are  solved  for 
the  special  case  of  constant  values  of  fluid  viscosity  and 
thermal  conductivity.  These  are  transformed  into  a  series 
of  ordinary  differential  equations  and  then  solved  with  the 
aid  of  the  Whirlwind  Digital  Computer  for  selected  entrance 
and  thermal -boundary  conditions . 

The  calculated  results  presented  in  the  main  text 
correspond  to  two  cases  at  an  entrance  Mach  number  of  2.8, 
with  and  without  heat  transfer.  They  include  recovery 
factor,  friction  coefficient,  mean  pressure-gradient  coeffi¬ 
cient,  heat -transfer  coefficient,  displacement  and  momentum 
thicknesses,  and  profiles  of  stagnation  temperature,  shear 
stress  and  heat-tran3fer  rate. 

The  results  for  tube  flow  are  compared  with  those  for 
plate  flow  with  corresponding  free-stream  and  thermal 
boundary  conditions.  Near  the  tube  inlet,  the  plate-flow 
solution  is  a  good  approximation  to  the  tube-flow  solution. 


♦This  investigation  is  sponsored  as  Contract  Number  N5- 
ori-07897  by  the  Office  of  Naval  Research  of  the  United  States 
Navy.  Availability  of  Digital  Computer  Laboratory  time  for  this 
problem  was  made  possible  also  by  the  Office  of  Naval  Research. 
Reproduction  of  this  report,  in  whole  or  in  part,  is  permitted 
for  any  purpose  by  the  United  States  Government. 

♦•Associate  Professor  of  Mechanical  Engineering. 

♦♦♦Professor  of  Mechanical  Engineering. 
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NOMENCLATURE 


a 

C 


D 

Fm 

£m 

fp 


h, i,  j,m,n 
hz 

Km 

k 


M 

P 

q 

R 

v,f,z 

r 

rd 

T 

t 

m 

u 

v 


fim 

8 

8* 

5 

V 

J 

e 


m 


inside  radius  of  tube 
constant,  see  eq.  (l) 

specific  heat  at  constant  pressure 

specific  heat  at  constant  volume 

inside  diameter  of  tube 

function  of  ?  ,  see  eqs .  ( 15 )  and  (l6) 

function  of  ,  see  eq.  (7) 

mean  pres sure -gradient  coefficient,  see  eq.  (A-19) 
local  true  friction  coefficient,  see  eq.  (A-16) 

any  positive  integer  for  exponent,  subscript,  or 
coefficient 

local  heat-transfer  coefficient,  see  eq.  (A-21) 

function  of  ,  see  eqs.  (14),  (17),  and  (18) 

ratio  of  specific  heats,  c  /c„ 

m  P  v 

coefficient  of  y  in  velocity  polynomial,  see 

eq.  (9) 

„  v  1/2 

Ma  ch  numbe  r ,  w/( kRT ) 
static  pressure 

rate  of  radial  heat  transfer  per  unit  circum¬ 
ferential  area 
perfect-gas  constant 
cylindrical  coordinates 
recovery  factor,  see  eq.  (A-l4) 
diameter  Reynolds  number,  w  D/V 

absolute  temperature 

coefficient  of  5m  in  temperature  polynomial, 
see  eq.  (ll) 

velocity  component  in  radial  direction 
velocity  component  in  axial  direction 
entrance  length 

coefficient  of  5  in  density  polynomial,  see 
eq.  (12) 

coefficient  of  in  pressure  polynomial,  see 
eq.  (10) 

boundary-layer  thickness 

displacement  thickness,  see  eq.  (A-l)  1/2 

modified  length  Reynolds  number,  2(z/DRrin) 

r  q  1  UyJ 

relative  distance  ratio,  ( r/&)‘~  J  A  5 
momentum  thickness,  see  eq.  (A-4) 
function  of  ,  see  eq.  (8) 
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thermal  conductivity 
absolute  viscosity 
kinematic  viscosity 
mass  density 

Prandtl  number,  c^jU/A 
shear  stress 
stream  function 


Subscripts : 


aw 

c 

0 

s 

T 

v 

w 

CD 

s 


re  f  e  rs 
refe  rs 
refers 
re  f  e  rs 
re  f  e  rs 
refe  rs 
refe  rs 
refers 
re  f  e  rs 


to  adiabatic-wall  conditions 

to  conditions  in  isentropic  core 

to  inlet  conditions 

to  stagnation  conditions 

to  temperature -boundary  layer 

to  velocity-boundary  layer 

to  conditions  at  tube  wall 

to  free-stream  conditions  for  plate  flow 

to  outer  edge  of  boundary  layer 


Primes  denote  differentiation  with  respect  to  J 


INTRODUCTION 


The  growing  importance  of  problems  involving  flow  of 
a  compressible  fluid  in  ducts  with  relatively  small  values 
of  the  length-to-diameter  ratio  has  been  recognized  in  many 
engineering  fields  during  the  past  decade.  In  these  short 
ducts,  the  flow  is  mainly  determined  by  the  nature  of  the 
development  of  the  boundary  layer  in  the  entrance  region, 
where  both  the  velocity  and  the  temperature  profiles  change 
rapidly  in  the  direction  of  flow. 

The  present  paper  presents  a  theoretical  solution  for 
the  case  of  a  compressible  fluid  flowing  at  either  subsonic 
or  supersonic  velocities  in  the  entrance  region  of  a  round 
tube,  when  a  laminar  boundary  layer  originates  at  the  tube 
inlet.  This  solution  is  obtained  from  the  integration  of 
the  partial  differential  equations  of  energy,  momentum,  and 
continuity,  for  selected  boundary  conditions,  for  a  gas 
following  the  perfect -gas  rule,  and  for  the  special  case  of 
constant  values  of  specific  heat,  viscosity,  and  thermal 
conductivity  of  air.  The  analysis  presented  here  is  not 
valid  when  the  laminar  boundary  layer  fills  the  entire  cross 
section  of  the  tube  or  when  the  boundary  layer  becomes 
turbulent . 

The  velocity  and  temperature  profiles  together  with  the 
variation  of  the  stream  properties  of  the  isentropic  core 
obtained  by  the  present  analysis  are  reported  by  Toong  and 

Kaye^*  .  rp^-pg  present  paper  presents  the  calculated  values 
of  the  recovery  factor,  friction  coefficient,  mean  pressure- 
gradient  coefficient,  and  heat-transfer  coefficient  for 
several  specific  boundary  conditions,  in  addition  to  the 
calculated  profiles  of  stagnation  temperature,  shear  stress, 
and  heat-transfer  rate  in  the  laminar  boundary  layer.  The 
variation  of  displacement  and  momentum  thicknesses  along 
the  length  of  the  tube  is  also  included. 


*Superscripts  refer  to  items  in  Bibliography . 


LITERATURE  SURVEY 

Both  analytical  and  experimental  studies  of  the  develop¬ 
ment  of  flow  in  the  entrance  region  of  a  roundtube  are  found 
in  the  literature.  Table  1  summarizes  briefly  the  analytical 
work  for  laminar  and  turbulent  flows  as  well  as  for  compressi¬ 
ble  and  incompressible  fluids*. 

Most  of  the  theoretical  analyses  are  limited  to  the  case 
of  a  laminar  flow  of  an  incompressible  fluid  in  a  round  tube . 

Hagenbach^,  Neumann-^,  and  Couette^  applied  merely  a  correction 
accounting  for  the  change  in  the  kinetic  energy  to  the  solution 

of  a  fully  developed  laminar  flow;  Schiller^,  and  Shapiro, 
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Siegel  and  Kline  used  the  integral  method;  Boussinesq  , 

Atkinson  and  Goldstein^,  and  Langhaar^  solved  the  partial 
differential  equations  of  momentum  and  continuity;  and 

Kays10  solved  the  energy  equation,  using  the  calculated  vel¬ 
ocity  profile  of  Langhaar. 

The  agreement  between  the  calculated  values  of  the  center- 
line  velocity  by  Boussinesq?  and  the  experimental  values  of 
Nikuradse11  is  within  3$  for  z/DRD0 >  0 .02 .  However,  this 

difference  is  increased  to  about  10$  at  smaller  values  of 
z/dRdq,  such  as  at  0.0075.  Atkinson  and  Goldstein^  used  the 
method  of  Boussinesq  to  extend  their  series  solution.  When 
these  two  solutions  are  joined  together  at  z/dRjjq  =  0 .000625, 

the  agreement  between  the  theoretical  center-line  velocity 
and  the  measured  value  of  Nikuradse  is  improved  to  within  3$ 
for  z/DR^  >  0.0075.  The  agreement  between  the  calculated 

values  of  the  velocity  profile  across  the  tube  section  and 
the  measured  values  of  Nikuradse  is  within  10$  for 
z/DRjjq  >  0.0075.** 

The  agreement  between  the  theoretical  value  of  the  center- 

Q 

line  velocity  of  Langhaar^  and  the  experimental  value  of 
Nikuradse11  is  within  10$  for  the  entire  entrance  length***; 
but  this  agreement  improves  to  better  than  3$  for  z/DRD0>0.01. 

*Those  analyses  dealing  with  diabatic  flows  with  fully  devel¬ 
oped  velocity  profiles  are  not  included  here . 

**See  Fig.  8l  on  p.  308  of  Golds tein°. 

***Entrance  length  obtained  by  the  integral  method  is  defined  as 
the  distance  at  which  cf/a=l,  and  that  obtained  from  the 
solution  of  the  boundary -layer  equations  is  defined  as  the 
distance  where  the  center-line  velocity  is  99$  of  the 
Poiseuille  value . 


For  small  values  of  z/DRd0  where  cf/a«l,  the  velocity 
profiles  assumed  in  the  integral  methods  of  S  chiller®  and 

Shapiro,  Siegel,  and  Kline®  agree  well  with  those  calculated 
by  the  more  exact  methods  of  solving  the  differential  equa¬ 
tions  of  momentum  and  continuity.  For  example,  for 

z/DR^q  <  0.000625,  the  velocity  profiles  of  Schiller®  agree 

with  those  of  Atkinson  and  Goldstein®  within  1 $  for  r/a  4  0.8 
and  within  5$  for  r/a  =0.9.  However,  for  such  small  values 

of  z/DR~^,  the  experimental  results  do  not  agree  with  the 
DO  ,  .  * 

theoretical  solutions . 

The  values  of  the  center-line  velocity  calculated  by 
Schiller  agree  within  10$  for  the  entire  entrance  length  with 
the  values  measured  by  Nikuradse .  The  agreement  gets  poorer, 
however,  as  z/DRqq  approaches  the  limiting  value  correspond¬ 
ing  to  the  entrance  length,  in  contrast  to  the  improved 
agreement  observed  between  the  velocities  calculated  from 
the  solution  of  the  boundary-layer  equations  and  the  measured 
values . 

The  entrance  lengths  calculated  by  the  integral  methods 

of  Schiller®  and  Shapiro  et  al .®  are  about  50$  smaller  than 
the  experimental  value  (0.043  DRD0)  of  Nikuradseli .  The 
entrance  lengths  obtained  by  Boussinesq' ,  Atkinson  and. 
Goldstein^,  and  Langhaar®  are  higher  than  the  Nikuradse ' s 
value  by  50$,  26$,  and  53$,  respectively. 

The  mean  pressure-gradient  coefficient  at  a  distance 
from  the  tube  inlet  greater  than  the  entrance  length  is 
calculated  from  an  equation  of  the  following  form: 

2  (p0  -  p)  //°0  wo  =  C  +  64  z/DRD0  (1) 

The  values  of  C  in  eq.  (l)  obtained  by  Boussinesq' ,  Atkinson 

and  Goldstein®,  and  Langhaar®  from  the  solutions  of  the 
boundary-layer  equations  agree  within  7$  with  the  mean 

experimental  value  (1.32)  of  Schiller^,  values  calcu¬ 

lated  by  the  integral  methods  of  Schiller®  and  Shapiro  et  al ,® 
and  by  the  method  of  kinetic -energy  correction  are  lower  by 
12$,  18$,  and  24$,  respectively,  than  the  experimental  value 

of  Schiller^  .  However,  it  is  important  to  note  also  that 
the  difference  between  the  mean  experimental  values  of 

Schiller^  and  those  of  Riemann^®  is  about  6$. 


A  comparison  of  the  theoretical  and  experimental  values 
of  the  mean  pressure -gradient  coefficient  at  a  distance  from 
the  tube  inlet  smaller  than  the  entrance  length  may  be  found 

in  the  work  of  Shapiro  et  al^  Over  the  range  of  the  experi¬ 
mental  values  of  10"’-'  <  z/DR^  <  10“^,  obtained  by  Kline  and 
14  ■LK~J  e; 

Shapiro  ,  the  calculated  values  of  Schiller^  and  Shapiro  et  al 

agree  within  60  and  20,  respectively,  of  the  experimental 
values.  For  z/DRdq > 0 .0075,  a  maximum  difference  of  70  is 

observed  between  the  calculated  values  of  Atkinson  and  Gold- 

Q  C 

stein0  and  Shapiro  et  al  .  As  z/DR  increases  toward  its 

DO 

limiting  value  corresponding  to  the  entrance  length,  this 
difference  is  reduced.  Fig.  5  of  reference  (6)  shows  that,  for 
z/DRjtjq  >  0  .0075,  the  mean  pressure -gradient  coefficients  calcu¬ 
lated  by  Atkinson  and  Goldstein  are  always  greater  than  those 


calculated  by  Shapiro  et  al 
lying  between  them. 


with  those  calculated  by  Langhaar-7 


As  z/DR__  approaches  zero,  the  experiments  of  Kline  and 

I  Ji 

Shapiro  show  that  the  ratio  of  the  mean  pressure -gradient 
coefficient  to  z/DR^  approaches  an  asymptotic  value  cf  1J.  .74. 

The  corresponding  asymptotic  values  obtained  from  the  analy¬ 
ses  of  Atkinson  and  Goldstein^,  Shapiro  et  al.^  and  Schiller^ 
are  0.20,  20,  and  60  higher,  respectively,  than  the  experimental 

value.  The  solution  of  Langhaar^  does  not  seem  to  have  the 
similar  asymototic  behavior  near  the  tube  inlet  as  others  . 

For  the  laminar  flow  of  a  compressible  fluid  in  the 
entrance  region  of  a  circular  tube ,  only  the  theoretical  analy¬ 
ses  of  Toong^-5  and  Toong  and  Shapiro^  are  available  .  Tocng^-5 
solved  the  partial  differential  equations  of  continuity, 
momentum,  and  energy,  for  selected  boundary  conditions  and  for 
the  special  case  of  constant  values  of  specific  heat,  viscosity, 
and  thermal  conductivity  of  air.  Toong  and  Shapiro^-0  used  an 
integral  method,  with  the  values  of  the  viscosity  and  thermal 
conductivity  of  air,  which  either  varied  according  to  the 
Sutherland's  equation  or  were  constant. 


As  for  the  case  of  an  incompressible  fluid,  the  integral 
method  gives  valid  results  only  when  .the  boundary-layer  thickness 
is  small  as  compared  with  the  tube  radius .  In  the  region  where 
6/ a  <0.3,  the  agreement  between  the  calculated  values  of  the 
mean  pressure -gradient  coefficient  by  the  two  analyses  quoted 
above  is  within  50 • 

17 

The  theoretical  analysis  of  Latzko  is  available  in  the 
study  of  the  turbulent  flow  of  an  Incompressible  fluid  in  the 
entrance  region  of  a  round  tube.  In  this  analysis,  the 


5 


developing  turbulent  boundary  layer  13  assumed  to  originate 
at  the  tube  inlet,  with  the  thickness  of  the  velocity¬ 
boundary  layer  equal  to  that  of  the  temperature -boundary 
layer. 


PROBLEM 

The  flow  model  of  the  present  paper  is  shown  in  Pig.  1. 
A  compressible  fluid  flows  steadily  in  the  entrance  region 
of  a  round  tube .  The  velocity  and  temperature  profiles  at 
the  tube  inlet,  or  station  0,  are  assumed  to  be  uniform.  A 
laminar  boundary  layer  is  assumed  to  grow  in  thickness  from 
the  inlet  station  until  at  some  point  downstream  the  laminar 
boundary  layer  either  fills  the  entire  cross  section  of  the 
tube  or  becomes  turbulent.  The  analysis  presented  in  this 
paper  is  not  valid  in  the  region  downstream  of  such  a  point. 

The  following  assumptions  are  used  in  the  analysis: 

1.  The  fluid  flow  inside  the  core  is  one-dimen¬ 
sional  and  isentropic. 

2.  Air  is  a  perfect  gas  with  constant  specific 
heats  . 

J> .  Negligible  external  body  forces  exist. 

4.  The  flow  is  rotationally  symmetric  about  the 
tube  axis . 

5.  The  usual  assumptions  for  boundary -layer 
problems  are  also  made;  viz.. 


u 

w 


«  1, 


<$V 


« 


1, 


and  since  cpya/A=0(l)  for  most  gases. 


<Tt 


~Z~ 


« 


1 


ANALYSIS 

The  analysis  of  the  present  problem  is  given  in  the 

earlier  paper  of  Toong  and  Kaye1.  Only  a  brief  outline  will 
be  given  below.  A  summary  of  some  important  physical  quanti¬ 
ties  is  listed  in  Appendix  A. 
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By  the  use  of  the  above  assumptions,  the  basic  equations 
of  continuity,  momentum,  and  energy  are  first  transformed  into 
the  following  equations  for  the  boundary-layer  flow: 


a  (ry>u)  +  ^5.  (r J>  „)  =  0 


-3E  =  o 

3  r 


/=<uf*  +  *£>-??*?  (4) 

AP  <»  H  +  -  H>  =  *  H  +  7  l  <  a  r|?)+/,(||f  (5) 

Next,  the  partial  differential  equations  (2),  (3),  (4), 
and  (5)  are  transformed  into  a  series  of  ordinary  differential 
equations  by  means  of  the  following  assumptions: 

1.  A  stream  function  %  ,  satisfying  the  continuity 
equation  (2)  such  that 


= 


H=-y°u  r'  J 

can  be  expressed  in  terms  of  the  following  convergent  series: 

$  =  -  “2  A,  “o  Z  5%  <?>  (7) 

where  J  >  the  modified  length  Reynolds  number,  and  >7  ,  the 
relative  distance  ratio,  are  the  two  new  independent  variables, 
and  are  defined  in  Nomenclature . 

2.  The  temperature  in  the  boundary  layer  can  be 
expressed  in  terms  of  the  following  convergent  series : 

T  =  T0 15°'\  (?)  (8) 
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5.  The  velocity  and  the  properties  of  the  fluid  in  the 
isentropic  core  are  expressed  in  terms  of  the  following 
series,  whose  coefficients  may  be  determined  from  the  condi¬ 
tions  of  steady  flow!5; 


iu 

WC  =  wo  +  2.5  km> 
p  =  p0  a  +  Z3m/g 

Tc  =  T0  (1  +  Z5mtm) 


(10) 


(11) 


A  =  A  U+Z5«m> 


(IS) 


For  the  special  case  of  constant  values  of  viscosity  and 
thermal  conductivity  of  the  compressible  fluid,  the  final 
transformed  equations  are  as  follows: 


F1  +  fi  Pi  =  0 


(1  Ar)  Q'l  +  <*[  +  (k-1)  V?  (f|)2  /  4  =  o 


(13) 


f"  +  f  f'  -  (m-1)  f '  F  +  mF '  f  =  K_,(?) 

ml®  1  m  1  m  din"5 

(l/c r)  0M  +  f  o'  -  (m-1)  f'  e  +m9'  f  4- (k-1 )  M2  f'  f'  /2Hl4) 

mim  1  m  m  0im 


=  K2m-2  (?] 


where 


V  8if  i 


(15) 


'  '  _ 

F  =  9  f  +  f.  9  ,  for  m  ^ 
m  l  m  1  m 


(16) 


8 


and  Kg  (  ? )  and  K 
for  example. 


2  m- 2 


(  >7 )  are  known  functions  of  >7  ; 


Kx  (?)  =  4  /51  /  kMQ2  +  4  F-l  "  f  \  (17) 

K2  (^)=  4  01  /cr  +(k-l)  M02A  (Fi.)2/2  "  (k-1)/^  FXA 


-4^fx  ox 


(18) 


Boundary  Conditions 


At  the  outer  edge  of  the  boundary  layer,  provided  that  an 
isentropic  core  exists,  it  may  be  shown  that1-^ 


m-1 


(f  '  )  =  2  k.  *  t  n 

m  c  z —  i  m-i-JL 


(19) 


and 


i=0 


(e  )  =  , 

'  m  e  m-1 


(20) 


with  kQ  =  o(q  =  t  =  1,  m  —  1,  2,  .... 

The  condition  that  the  velocity  vanish  at  the  tube  wall 
requires  that 


t  (0)  =  0 

m 


(21) 


f  (0)  =  0 
m 

The  condition  of  an  adiabatic  wall  requires  that 


(22) 


®m  <°)  ®  0 


(25) 


and  the  condition  of  a  wall  maintained  at  an  arbitrary  tempera¬ 
ture  distribution  requires  that 

m-1 

Z  5  e„  (0) = 


(24) 
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RESULTS 

The  transformed  boundary -layer  equations  (15)  and  ( 14 ) 
for  the  special  case  of  constant  values  of  viscosity  and 
thermal  conductivity  of  air  were  solved  first  with  the  aid 
of  the  M.I.T.  Differential  Analyzer  and  more  recently  with 
the  aid  of  the  Whirlwind  I  Digital  Computer  for  better 
accuracy.  The  solutions  obtained  from  the  Whirlwind  Compu¬ 
ter  agree  with  those  obtained  from  the  Differential  Analyzer 
to  the  last  significant  figure  of  the  latter  solutions .  The 
results  reported  in  this  paper  are  obtained  from  the  Whirl¬ 
wind  Computer  and  are  tabulated  in  Appendix  B. 

Several  solutions  of  the  equations  were  obtained  for 
different  sets  of  entrance  and  thermal -boundary  conditions, 
but  only  two  cases  are  presented  in  this  section.  Both  of 
these  cases  correspond  to  an  entrance  Mach  number  of  2.8; 
case  A  is  for  adiabatic  flow  and  case  C  is  for  flow  with 

heat  transfer  to  the  compressible  fluid  at  0.  =  3 .  The 

W 

velocity  and  temperature  profiles  together  with  the  variation 
of  the  stream  properties  in  the  isentropic  core  corresponding 
to  these  two  cases  have  been  reported  by  Toong  and  Kaye* . 
Calculated  values  of  recovery  factor,  friction  coefficient, 
mean  pressure -gradient  coefficient,  heat-transfer  coefficient, 
displacement  and  momentum  thicknesses,  and  profiles  of  stag¬ 
nation  temperature,  shear  stress,  and  heat-transfer  rate  are 
given  here.  These  physical  quantities  can  be  computed  from 
the  equations  given  in  Appendix  A. 

Boundary -Layer  Thickness -Ratios 

The  boundary-layer  thickness  or  <£,p  is  related  through 

eq.  (A-3)  in  Appendix  A  to  the  corresponding  value  of  the 
relative  distance  ratio  at  which  the  velocity  or  tempera¬ 
ture  inside  the  boundary  layer  reaches  a  certain  arbitrary 
percentage  of  the  core  velocity  or  temperature,  respectively. 
The  upper  part  of  Pig.  2  shows  the  variation  of  the  thickness 
ratios  for  both  the  velocity-and  temperature -boundary  layers, 
when  the  first  terms  of  the  series  solutions  for  the  velocity 
and  temperature  in  the  boundary  layer  are  equal  to  99$  of  the 
core  velocity  and  core  temperature,  respectively. 

Pig.  2  Indicates  that  the  thickness  of  the  velocity¬ 
boundary  layer  is  smaller  than  that  of  the  temperature - 
boundary  layer,  for  given  values  of  the  entrance  Mach  number 
and  given  thermal  conditions  at  the  tube  wall.  This  result 
is  also  predicted  from  a  study  of  the  basic  equations  of 
motion;  i.e.,  the  ratio  of  the  thickness  of  the  velocity¬ 
boundary  layer  to  that  of  the  temperature -boundary  layer  is 
shown  to  be  of  the  order  of  the  square  root  of  the  Prandtl 
number  of  the  fluid  flowing  inside  the  tube . 
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For  given  values  of  the  entrance  Mach  number,  the  effect 
of  heating  the  fluid  at  the  tube  wall  is  to  increase  the 
boundary -layer  thicknesses . 

Shown  for  comparison  in  the  upper  part  of  Fig.  2  are  also 
the  boundary-layer  thickness-ratios  predicted  for  an  adiabatic 
laminar  flow  of  a  compressible  fluid  over  a  flat  plate,  with 
the  free -stream  conditions  identical  with  the  entrance  condi¬ 
tions  of  the  corresponding  tube  flow.  Near  the  tube  inlet,  it 
is  observed  that  the  plate-flow  solution  is  a  good  approxi¬ 
mation  to  the  tube-flow  solution. 

The  middle  part  of  Fig.  2  shows  the  variation  of  the 
displacement  and  momentum  thicknesses,  defined  by  eqs .  (A-l) 

and  (A-4),  for  tube  flows  at  Mq-2.8  with  and  without  heat 

transfer.  Two  dashed  curves  are  also  included  to  show  the 
variation  of  these  thicknesses  for  an  adiabatic  laminar  flow 
over  a  flat  plate  at  a  free-stream  Mach  number  of  2.8.  Near 
the  tube  inlet,  the  curves  for  tube  flow  are  almost  linear 
with  respect  to  the  modified  length  Reynolds  number,  thus 
demonstrating  again  the  characteristics  of  a  plate  flow. 

The  ratio  of  displacement  thickness  to  momentum  thick¬ 
ness  is  shown  in  the  lower  part  of  Fig.  2.  For  a  given 
entrance  Mach  number  and  thermal  condition  at  the  tube  wall, 
this  ratio  is  almost  independent  of  the  modified  length 
Reynolds  number.  Shown  for  comparison  are  also  two  curves 

for  laminar  flow  over  a  flat  plate  at  Mt0=  2  .8,  with  and 

without  heat  transfer.  For  plate  flows,  the  ratio  of  these 
thicknesses  remains  constant. 

Region  of  Validity  of  Present  Analysis 

Since  the  present  analysis  is  valid  only  when  the  boundary 
layer  does  not  fill  the  entire  cross  section  of  the  tube,  the 
calculated  results  should  not  be  extended  beyond  a  certain 
maximum  value  of  the  modified  length  Reynolds  number.  Since 
the  present  analysis  is  valid  only  when  the  boundary  layer 
remains  laminar,  the  calculated  results  should  not  be 
extended  beyond  a  certain  value  of  the  transitional  length 
Reynolds  number.  Fig.  5  is  constructed,  on  the  basis  of  these 
limitations,  to  show  the  region  of  validity  of  the  present 
analysis  for  an  entrance  Mach  number  of  2.8.  Two  values  are 
selected  for  both  the  maximum  modified  length  Reynolds  number 
and  the  transitional  length  Reynolds  number.  The  values  of 
the  length-to-diameter  ratios  below  the  straight  line  of  a 
given  modified  length  Reynolds  number  and  the  rectangular 
hyperbola  of  a  given  transitional  length  Reynolds  number  are 
in  the  valid  region  of  the  present  analysis.  It  is.  observed 
that  the  limiting  condition  at  low  entrance  diameter  Reynolds 
numbers  is  the  value  of  the  modified  length  Reynolds  number. 
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while  the  limiting  condition  at  high  entrance  diameter 
Reynolds  numbers  is  the  value  of  the  transitional  length 
Reynolds  number. 

Profiles  of  Ts/Tsc,  T/Tw,  qD/*0  TQ,and  q/c^ 

Pig.  4  presents  for  the  adiabatic  flow  at  Mq=2.8 

the  profiles  of  the  stagnation  temperature,  shear  stress, 
and  heat-transfer  rate  across  the  laminar  boundary  layer 
at  three  different  axial  positions.  They  are  calculated 
from  eqs  .  (A -6),  (A-ll),  and  (A-12),  respectively.  The 

first  two  profiles  for  Tg/Tgc  and  T/<2W  are  shown  at 
5=0,  0.03,  and  0.05,  and  the  third  profile  for  qD/^0T0 
is  shown  at  5  “0.01,  0.03,  and  0.05. 

Pig.  5  presents  the  profiles  of  the  stagnation 
temperature,  shear  stress,  and  heat-transfer  rate  for  flow 

with  heat  transfer  to  the  compressible  fluid  at  0y=  3  and 
Mq=2.8.  They  are  calculated  from  eqs.  (A-6),  (A-ll),  and 
(A-13),  respectively.  Note  the  difference  in  the  ordinate 
representing  the  heat-transfer  rate  in  Pigs.  4  and  5.  Also 
note  that  in  Fig.  5,  all  the  three  profiles  are  shown  at 
5  =  0,  0.03,  and  0.05. 

Near  the  outer  edge  of  the  boundary  layer  adjacent  to 
the  core,  a  slight  bump  is  found  at  the  larger  values  of  5 
for  the  profiles  shown  in  Pigs.  4  and  5.  Similar  bump3  are 
observed  in  the  velocity  and  temperature  profiles  in  the 

earlier  paper  of  Toong  and  Kaye^ .  These  bumps  have  actually 
been  measured  for  supersonic  flow  of  air  in  a  round  tube 
by  Kaye,  et  al.1®'^ 

The  bump  in  the  stagnation-temperature  profile  of 
Fig.  4  can  be  checked  by  evaluating  the  integral  of  the 
enthalpy  flux  with  respect  to  radius  for  various  sections 
at  $  —  0,  0.03,  and  0.05.  These  fluxes  were  found  to  remain 
constant  within  1$. 

It  is  to  be  noted  that  the  profiles  at  5  —  0  not  only 
show  the  solution  at  the  tube  inlet  but  also  represent  the 
solution  of  the  corresponding  plate  flow.  In  fact,  for 
adiabatic  flow  at  Mq  =  2.8,  the  agreement  between  the  tube 

solution  at  5  =0  and  the  solution  of  Emmons  and  Brainerd20 
is  within  1$.  Since  the  solution  at  5=0  is  a  first 
approximation  to  the  series  solution  of  the  tube  flow  at 
some  finite  distance  from  the  inlet,  it  is  demonstrated 
once  again  that  near  the  tube  inlet,  the  tube  flow  behaves 
like  a  plate  flow  with  corresponding  free-stream  conditions. 
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Recovery  Factor,  Friction  Coefficient,  Mean  Pressure- 
Gradient  Coefficient,  and  Heat-Transfer  Coefficient 

Fig*  6  shows  the  variation  of  recovery  factor,  friction 
coefficient,  mean  pres sure -gradient  coefficient,  and  heat- 
transfer  coefficient  with  respect  to  the  modified  length 
Reynolds  number. 

The  recovery  factor,  defined  by  eq.  (A-l4)  and  calculated 
from  eq.  (A-15)  is  plotted  at  the  top  of  Fig.  6  for  the  adia¬ 
batic  flow  at  Mq=2.8.  The  recovery  factor  increases  almost 

linearly  from  0.846  to  0.855  as  the  value  of  _5  increases  from 
zero  to  0.04. 

The  theoretical  analysis  of  Emmons  and  Brainerd^0  for  an 
adiabatic  laminar  flow  of  a  compressible  fluid  over  a  flat 
plate  indicates  that,  as  a  good  first  approximation,  the 
recovery  factor  is  equal  to  the  square  root  of  the  Prandtl 
number  of  the  compressible  fluid.  In  Fig.  6,  the  line  repre¬ 
senting  the  square  root  of  the  Prandtl  number  (0.74)  used  in 
the  present  calculation  is  also  shown  for  comparison.  Thus 
the  tube -flow  solution  near  the  inlet  is  found  to  be  one  per 
cent  less  than  the  approximate  plate -solut ion . 

The  local  true  friction  coefficient,  defined  by  eqs . 

(A-l6)  and  (A-17)  and  calculated  from  eq.  (A-l8),  is  shown 
in  Fig.  6  for  flows  at  Mq  =  2.8,  with  and  without  heat  trans¬ 
fer.  The  ordinate  is  in  the  form  of  the  product  of  the 
friction  coefficient  and  the  square  root  of  the  length 
Reynolds  number.  For  laminar  flow  over  a  flat  plate,  this 
product  is  constant;  while  for  laminar  flow  in  a  round  tube, 
this  product  decreases  almost  linearly  by  about  15  to  20  $ 
as  the  value  of  5  increases  from  zero  to  0.04.  Shown  for 
comparison  in  Fig.  6  is  also  the  curves  for  laminar  flow  over 
a  flat  plate  at  a  free-stream  Mach  number  of  2.8,  with  and 
without  heat  transfer.  Again,  it  is  observed  that  near  the 
tube  Inlet,  the  plate -flow  solution  is  a  good  approximation 
of  the  tube-flow  solution. 

Fig.  6  also  Indicates  that  the  effect  of  heating  is  to 
decrease  the  friction  coefficient  for  the  same  entrance 
diameter  Reynolds  number  and  at  the  same  length  Reynolds 
number. 

The  mean  pressure -gradient  coefficient,  defined  according 
to  eq.  (A-19)  and  calculated  from  eq.  (A-20;  is  shown  in  Fig.  6 
for  flows  at  MQ  =  2.8,  with  and  without  heat  transfer.  Here 

it  is  observed  that  the  effect  of  heating  the  compressible 
fluid  at  the  tube  wall  on  the  mean  pressure -gradient  coeffi¬ 
cient  is  just  opposite  to  that  on  the  local  true  friction 
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coefficient;  viz.,  heating  gives  higher  values  of  the 
pressure -gradient  coefficient  but  lower  values  of  the  friction 
coefficient.  This  apparent  discrepancy  is  due  to  the  differ¬ 
ent  ways  the  two  coefficients  are  defined.  The  pressure- 
gradient  coefficient  is  related  to  the  change  in  static 
pressure  along  the  length  of  the  tube.  Using  the  knowledge 
developed  in  the  one -dimensional  flow  of  a  compressible 
fluid  inside  a  duct  of  constant  area  that  the  effects  on 
pressure  change  due  to  heating  and  friction  are  in  the  same 
direction,  it  seems  justified  that  the  value  of  the  pressure  - 
gradient  coefficient  is  increased  when  the  fluid  is  heated 
at  the  tube  wall.  On  the  other  hand,  the  friction  coefficient 
is  related  to  the  radial  velocity  gradient  at  the  tube  wall . 
When  the  fluid  is  heated,  the  boundary -layer  thickness  is 
increased,  the  radial  velocity  gradient  at  the  wall  is 
decreased,  and  therefore  the  friction  coefficient  is 
decreased . 

At  the  bottom  of  Fig.  6  is  shown  a  plot  of  the  product 
of  the  Stanton  number  and  the  square  root  of  the  length 
Reynolds  number,  for  laminar  flow  with  heat  transfer  at 

0=3  and  =  2.8.  The  heat-transfer  coefficient  is 
w  ^  u 

defined  by  eq.  (A-21),  and  the  Stanton  number  is  calculated 
from  eq.  (A-22).  It  is  observed  that  this  product  Is 
almost  independent  of  the  modified  length  Reynolds  number. 

In  the  same  plot  Is  shown  also  a  line  for  the  laminar 
flow  of  a  compressible  fluid  over  a  flat  plate  at  1^0,=  2.8 

and  Tw/T<d=3.  The  two  solutions  agree  within  2$. 


/ 
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OP  THEORETICAL  ANALYSES  POP,  PLOW  IN  THE  ENTRANCE  REGION  OF  A  ROUND  TUBE 
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APPENDIX  A 


Summary  of  Important  Physical  Quantities 

1.  Displacement  Thickness,  6* 

r°- 

fovc  n  <a  -  ■5*)  =  Awo  n  &2  -  271  /  </owo  -/•»)  rdr  (A-l) 

or  2(£L)  -  =  H  5 /  (l  -  j*.)  (A-2) 

where  ^  =  [l  -  (1  -  _4-)2  ]  /  *5  (A-3) 

2.  Momentum  Thickness,  ^ 

7t  (a  -^)8  s  Tl  a2  -  2/t  /  (w0  -  w)_/>w  rdr  (A-4) 

2  A  y“"'5 

(-B-)  -  <^)  =  4  5/  (1  -  5  )  d?  (A-5) 


or  2 


3.  St  agnation -Temperature  Profile,  T_/T 

s  sc 


Is-  =  =T(l+i-§-i  M2)  /  T0(l  +  i-^i  M2) 

sc  s  0 


with 


M2 

w 


2  ,  T. 


r©  ©  (t 


is 

w. 


w 

w. 


=  1  +  Z5m  >s, 


2  j"  +  n  '  2  e.  r 


m  n 


2(1  +  +  ZSlti)<1  +  Z5jkj) 


T  _  ^r-  _  m  -  1 

-  Z5  » 


m 


(A-6) 

(A-7) 

(A-8) 

(A-9) 


and 


(A-10) 


lb 


4.  Shear-Stress  Profile,  't /Zy 


T  _  /  c)v  )/fdv  \ 

Tw  (  2>v)/  U  r  yw 


r25m  +  "  '  2  (em  fn  +  ®i  fn) 


Z5 


m  +•  n  -  2 


(®m  *n  +  em  fn)w 


5.  Heat -Transfer-Rate  Profile,  q  dAqTq  or  q/qw 


A3^-  = 

Ao  To 


r 

a 


(2.5 


m  -  2 


©  ) 
m 


and 


Q  __ 


r  (  Z  5"  '  1  O 


<2  5 


m  -  1  q  '  \ 


m/w 


6.  Recovery  Factor,  r 


r  = 


or  r  = 


(T  -  T  )  /  (T  -  T  ) 
v  aw  c;  '  v  sc  c 1 


-  1  8m  (0)  -  (1  tZ5n  tn) 

(k-1)  M2/2  -  Z3n  tn 


7.  Local  True  Friction  Coefficient,  fF 


'  *  V  <s  'Vo) 


with 


=  >“( l?)w 


or  tz  = 


(Vv>1/2  Z  5m+n  '  S(®m  fn  +  «i  fA>v 


(A-ll) 


(A-12) 

(A-13) 

(A-14) 

(A-15) 


(A-16) 

(A-17) 


2 


(  1 + Z5m/3m) 


(A-18) 
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8.  Mean  Pressure-Gradient  Coefficient,  4  f  (z/d) 
4  fp  (z/D)  =  (p  -  P0)  /\  f0  v* 

or  4  fp  (Z/D)=S(2>5n/6n  -1)AM2 

9.  Local  Heat -T rans f e r  Coefficient,  h2 

h  EE  q/(T  -  T  ) 
z  '  v  aw 

h 

or  z 

Cp/°0W0 


(A-19) 

(A-20) 

(A-21) 


^  ($)1/8Z3m  '  1(e:>w/Z3"  '  1  [(eA  -  (8^  j  (A-K) 
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APPENDIX  B 

TABULATION  OF  RESULTS 

The  transformed  boundary-layer  equations  (13)  and  (14) 
for  the  special  case  of  constant  values  of  fluid  viscosity 
and  thermal  conductivity  were  solved  with  the  aid  of  the 
Whirlwind  I  Digital  Computer  for  the  six  cases  shown  below  in 
Table  B-l.  The  Prandtl  number  of  the  fluid  was  assumed  con¬ 
stant  at  0.74.  The  series  solution  of  each  of  these  cases  wa3 

carried  up  to  terms  containing  \  with  the  value  of  m 
indicated  in  the  table . 


TABLE  B-l 


Entrance  Thermal  Condition 


Case 

Mach  Number 

at  Tube  Wall 

m 

A 

2.8 

Adiabatic 

5 

B 

2.8 

Cooling,  9W=  2 

2 

C 

2.8 

Heating,  9W=  3 

3 

D 

2 

Adiabatic 

2 

E 

2 

Cooling,  9W= 1.445 

2 

F 

2 

Heating,  9w=2.1o7 

2 

The  functions  fn  ,  Fn  ,  0n  ,  and  their  derivatives  are 

tabulated  in  Table  B-2,  and  the  quantities  ,  tn  ,  (Xn,  and 
k  in  Table  B-3.  The  physical  quantities  listed  in  Appendix  A 

can  be  computed  from  these  tables  . 
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TABLE  B  -2 

Values  of  fn,  Fn,  Fn',  9n,  and  0^ 
(1)  Case  A,  M0=2.8,  Adiabatic 


V 

fl 

pl 

pl' 

el 

-V 

0.0 

.1 

.2 

.3 

.4 

0.00000 
.  208 
.  834 

.  1881 
.  3357 

0.00000 
.  9678 
.19353 
.29020 
.38668 

0.96779 

.96772 

.96725 

.96597 

.96348 

2.3269 

2.3242 

2.3160 

2.3025 

2.2835 

0.00000 
.  5433 

.10861 

.16269 

.21633 

.3 

.6 

.7 

.8 

.9 

.  5272 
.  7636 
.10466 
.13778 
.17592 

.48284 

.57849 

.67341 

•76733 

.85997 

.95937 

.95323 

.94468 

.93334 

.91885 

2.2592 

2.2297 

2.1951 

2.1556 

2.1115 

.26920 

.32087 

.37082 

.41847 

.46314 

1.0 

1.1 

1.2 

1.3 

1.4 

.21931 

.26820 

.32283 

.38355 

.45062 

.95099 

1.04003 

1.12671 

1.21063 

1.29140 

.90091 

.87926 

.85370 

.82411 

.79049 

2.0631 

2.0108 

1.95514 

1.89660 

1.83582 

.50413 

.54070 

•57213 

.59771 

.61681 

1.5 

1.6 

1.7 

1.8 
1.9 

.52435 

.60506 

.69306 

.78861 

.89197 

1.36860 

1.44186 

1 .51081 
1.57515 
1.63461 

.75292 

.71162 

.66694 

.61936 

.56948 

1.77347 

1.71028 

1 .64699 
1.58435 
1.52310 

.62893 

.63369 

.63091 

.62064 

.60312 

2.0 

2.1 

2.2 

2.3 

2.4 

1.00330 

1.12273 

1.25028 

1.38585 

1.52927 

1.68899 

1.73819 

1.78216 

1 .82097 
1.85476 

.51803 

.46582 

.41373 

.36266 

.31349 

1.46395 

1.40753 

1 .35440 
1.30502 

1 .25972 

.57885 

.54858 

.51322 

.47386 

.43170 

2.5 

2.6 

2.7 

2.8 
2.9 

1.68020 

1.83824 

2.0029 

2.1735 

2.3494 

1.88376 

1.90827 

1.92868 

1.94539 

1.95883 

.26703 

.22397 

.18484 

.15002 

.11966 

1 .21873 
1.18214 
1.14993 
1.12197 
1.09803 

.38797 

.34389 

.30059 

.25905 

.22011 

3.0 

3.1 

3.2 

3.3 

3.4 

2.5300 

2.7147 

2.9026 

3.0934 

3.2863 

1.96947 

1.97773 

1.98402 

1.98874 

1.99220 

.09376 
.  7213 
.  5447 
.  4036 
.  2934 

1.07784 

1 .06104 

1 .04727 
1.03615 

1 .02728 

.18435 
.15220 
.12385 
.09932 
.  7850 
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TABLE  B-2  (cont'd) 


(i)  Case  A,  Mq=2.8,  Adiabatic  (cont'd) 


V 

*1 

F1 

V 

el 

-V 

3.5 

3.4811 

1.99469 

0 .02092 

1.02055 

0.06115 

3.6 

3.6772 

1.99645 

.  1463 

1.01495 

.  4695 

3.7 

3.8744 

1.99767 

.  1003 

1.01085 

.  3550 

3.8 

4.0724 

1.99850 

.  0674 

1.00777 

.  2646 

3-9 

4.2709 

1.99906 

.  444 

1.00549 

.  1944 

4.0 

4.4700 

1.99942 

.  287 

1 .00585 

.  1407 

4.1 

4.6693 

1.99965 

.  182 

1.00265 

.  1005 

4.2 

4.8688 

1.99979 

•  113 

1.00179 

.  0705 

4.3 

5.0685 

1.99988 

.  069 

1.00120 

.  488 

4.4 

5.2683 

1.99993 

41 

1.00079 

.  333 

4.5 

5.4682 

1.99997 

24 

1.00052 

.  224 

4.6 

5.6681 

1.99998 

14 

1.00053 

.  148 

4.7 

5.8680 

1.99999 

.  08 

1.00021 

.  097 

4.8 

6.0680 

2.0000 

4 

1.00015 

.  62 

4.9 

6.2680 

2.0000 

2 

1.00008 

.  39 

5.0 

6.4679 

2.0000 

1 

1.00005 

.  25 

5.1 

6 .6679 

2.0000 

1 

1.00005 

.  15 

5.2 

6.8679 

2.0000 

0 

1.00002 

.  09 

5.3 

7.0679 

2.0000 

0 

1.00001 

5 

5.4 

7.2679 

2.0000 

0 

1.00001 

3 

5.5 

7.4679 

2.0000 

0 

1.00001 

2 

5.6 

7.6679 

2 .0000 

0 

1.00001 

1 

5.7 

7.8679 

2.0000 

0 

1.00000 

1 

5.8 

8.0679 

2.0000 

0 

1.00000 

0 

5.9 

8.2679 

2.0000 

0 

1.00000 

0 

6.0 

8.4679 

2.0000 

0 

1.00000 

0 

O  r-l  OJ  ITWO  t^OO  CTi 


I 


V 
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TABLE 

B-2  (cont'd) 

(i) 

Case  A,  Mq  = 

2.8,  Adiabatic  (cont'd) 

1 

i 

f2 

P2 

v 

92 

1 

o'  ■ 

W2 

0.0000 

0  .0000 

7  .806 

0.8647 

0.0000 

.  166 

.8206 

8.605 

.8803 

.2893 

.  686 

1 .7205 

9.390 

.9182 

.4475 

.1597 

2.698 

10.148 

.9655 

.4776 

.2937 

3.749 

10 .864 

1.0096 

.3853  ; 

I 

.4745 

4  .869 

11.522 

1.0387 

.1795  | 

.7066 

6.051 

12.102 

1.0421 

-  .1272  I 

.9946 

7.286 

12.589 

1.0105 

-  .5179  ! 

1.3440 

8.565 

12.962 

0.9364 

-  .9720  ! 

1.7603 

9.874 

13.205 

.8148 

-1.4647  j 

1.0 

2.250 

11.201 

13.301 

.6431 

-1.9676 

1.1 

2.820 

12.529 

13.234 

.4220 

-2.450 

1.2 

3.477 

13-842 

12.994 

.1549 

-2.879 

1.3 

4.229 

15.122 

12.574 

-  .1511 

-3.224 

1.4 

5.084 

16.350 

11.971 

-  .4860 

-3.453 

1.5 

6.050 

17.510 

11 .191 

-  .8370 

-3.543 

1 .6 

7.134 

18.583 

10.244 

-1.1892 

-3.473 

1.7 

8.342 

19.554 

9.150 

-1.5260 

-3.234 

1.8 

9.680 

20.41 

7.937 

-1.8304 

-2 .827 

1.9 

11.149 

21.14 

6.639 

-2.086 

-2.263 

2.0 

12.749 

21.74 

5.297 

-2.278 

-1.5641 

2.1 

14.475 

22  .20 

3.956 

-2.396 

-0.7637 

2.2 

16.520 

22.53 

2  .663 

-2.429 

.0976 

2.3 

18.269 

22.73 

1.4630 

.-2.376 

.9742 

2.4 

20.31 

22.83 

0.3952 

-2.235 

1 .8194 

2.5 

22.41 

22.82 

-0 .5084 

-2.014 

2.590 

2.6 

24.55 

22.73 

-1.2271 

-1.7212 

3.247 

2.7 

26.71 

22.58 

=1.7525 

-1.3693 

3.765 

2.8 

28.87 

22.39 

-2  .089 

-0.9733 

4.127 

2.9 

31.00 

22.17 

-2 .252 

-  .5492 

4.328 

3.0 

33-08 

21.94 

-2.267 

-  .1129 

4.374 

3.1 

35.11 

21.72 

-2.165 

.3209 

4.280 

3.2 

37.06 

21.51 

-1.9788 

.7393 

4.071 

3-3 

38.95 

21.32 

-1.7405 

1.1320 

3.770 

3.4 

40  .76 

21.16 

-1.4788 

1.4914 

3 .4o8 
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TABLE  B-2  (contd.) 


(i)  Case  A,  MQ  = 


?  f2  P2 


3.5 

42.50 

21.03 

3.6 

44.19 

20.92 

3.7 

45.82 

20.83 

3.8 

47.40 

20.77 

3-9 

48.94 

20.72 

4.0 

50.46 

20.68 

4.1 

51.94 

20.66 

4.2 

53.41 

20.64 

4.3 

54.86 

20.63 

4.4 

56.30 

20.62 

4.5 

57.73 

20.61 

4 .6 

59.16 

20.61 

4.7 

60.58 

20.61 

4.8 

62.00 

20.61 

4.9 

65.42 

20.61 

5.0 

64.83 

20.61 

5.1 

66.25 

20.61 

5.2 

67.66 

20.61 

5.3 

69.07 

20.61 

5.4 

70.49 

20.61 

5.5 

71.90 

20.61 

5.6 

73.31 

20.61 

5.7 

74.73 

20.61 

5.8 

76.14 

20.61 

5-9 

77.55 

20.61 

6.0 

78.96 

20.61 

2.3,  Adiabatic  (contd.) 


0  1 
y2 


-1 .2170 

1.8125 

3.010 

-0.9718 

2.093 

2.601 

-  .7541 

2.333 

2.201 

-  .5694 

2.534 

1.8246 

-  .4186 

2.699 

1.4831 

-  .2999 

2.832 

1.1825 

-  .2095 

2.937 

0-9253 

-  .1428 

3.019 

.7108 

-  .0949 

3.081 

.5363 

-  .  616 

3.127 

.3974 

-  .  390 

3.161 

.2894 

-  .  241 

3.186 

.2071 

-  .  145 

3.203 

.1458 

-  .  085 

3.215 

.1009 

-  .  48 

3.224 

.0687 

-  .  26 

3.230 

.  460 

-  .  13 

3.233 

.  303 

-  .  06 

3.236 

.  197 

-  .  2 

3.237 

.  127 

0 

3.238 

.  080 

1 

3.239 

.  51 

2 

3.239 

.  32 

2 

3.240 

.  20 

2 

3.240 

•  13 

2 

3.240 

.  09 

3  5.240  .  6 


O  1-1  CVJ  KV4-  mvo  t>-co  o\ 


25 


TABLE  B-2  (cont'd) 


(i)  Case  A,  Mo=2.8>  Adiabatic  (cont'd) 


V 


1.0 

1.1 

1.2 

1.3 

1.4 

1.5 

1.6 

1.7 

1.8 

1.9 


f3 

P3 

V 

03 

V 

0.000 

0.000 

-22.0 

-0.918 

0.00 

-.  39 

-1.643 

-10.86 

-  .784 

2.68 

-.123 

-2.16 

0.506 

-  .388 

5.20 

-  .207 

-1.534 

12.10 

.243 

7.32 

-  .243 

0.263 

23.9 

1.055 

8.82 

-.183 

3.24 

35.7 

1.979 

9.52 

.024 

7.41 

47.6 

2.93 

9.32 

.432 

12.75 

59.2 

3.81 

8.18 

1.099 

19.25 

70.6 

4.54 

6.16 

2.09 

26.9 

81.4 

5.02 

3.43 

3-47 

35.5 

91.4 

5.21 

0.251 

5.32 

45.1 

100.2 

5.07 

-3-02 

7.72 

55.5 

107.6 

4.61 

-5.94 

10.76 

66.5 

112.9 

3.91 

-8.06 

14.51 

78.0 

115.7 

3.05 

-8.90 

19-04 

89.6 

115.3 

2.18 

-8.10 

24.4 

100.9 

111.1 

1.490 

-5-38 

30.7 

111  .6 

102.4 

1.171 

-0.675 

37.8 

121.2 

88.7 

1.416 

5.87 

45.8 

129.2 

69.5 

2.39 

13.85 

2.0 

54.4 

135.0 

45.2 

4.21 

22.7 

2.1 

63.5 

138.1 

16.53 

6.93 

31.5 

2.2 

72.7 

138.2 

-15.03 

10.49 

39.5 

2.3 

81.7 

135.1 

-47.2 

14.77 

45.7 

2.4 

90.1 

128.8 

-77.4 

19.55 

49.5 

2.5 

97.5 

119.8 

-102.6 

24.6 

50.3 

2.6 

103.4 

108.5 

-120.6 

29.5 

47.9 

2.7 

107.7 

95.9 

-129.9 

34.0 

42.4 

2.8 

110.2 

82.9 

-130.1 

37.9 

34.2 

2.9 

110.8 

70.2 

-122.1 

40.8 

24.0 

3.0 

109.7 

58.6 

-108.1 

42.7 

12.67 

3.1 

107.3 

48.7 

-  90.2 

43.3 

1.090 

3.2 

103.8 

40.6 

‘  71.1 

42.9 

-9.92 

3.3 

99.7 

34.4 

-  52.8 

41.4 

-19.65 

3.4 

95.3 

30.0 

-  36.8 

39.0 

-27.6 
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TABLE  B-2  (cont'd) 


(i)  Case  A,  MQ  =  2.8,  Adiabatic  (cont'd) 


V 

f3 

p? 

P3 

1 

95 

Q  1 

y3 

3.5 

90.9 

27.0 

-23 

.8 

36.0 

-33.4 

3.6 

86.8 

25.1 

-14 

.03 

32.4 

-37.1 

3.7 

83.3 

24.1 

-  7 

.24 

28.6 

-38.6 

3.8 

80.5 

23.6 

-  2 

.90 

24.7 

-38.4 

3.9 

78.3 

23.4 

-  0 

.400 

21.0 

-36 .6 

4.0 

76.8 

23.5 

.817 

17.47 

-33.7 

4.1 

76.0 

23.6 

1 

.232 

14.27 

-30.1 

4.2 

75.8 

23.7 

1 

.202 

11 .46 

-26.1 

4.3 

76.1 

23.8 

0 

.971 

9.06 

-22  .0 

4.4 

76.9 

23.9 

.687 

7.05 

-18.16 

4.5 

78.1 

23.9 

.430 

5.42 

-14.62 

4.6 

79-5 

24.0 

.231 

4.11 

-11.50 

4.7 

81.2 

24.0 

•  095 

3.10 

-  8.85 

4.8 

83.1 

24.0 

.  11 

2.33 

-  6.67 

4.9 

85.1 

24.0 

- 

.  32 

1.751 

-  4.93 

5.0 

87.2 

24.0 

_ 

.  49 

1.329 

-  3.56 

5.1 

89.3 

24.0 

- 

.  51 

1.027 

-  2.53 

5.2 

91.5 

24.0 

- 

.  45 

0.814 

-  1.762 

5.3 

93.8 

24.0 

- 

.  36 

.668 

-  1.204 

5.4 

96.1 

24.0 

- 

.  27 

.568 

-  0.807 

5.5 

98.4 

24.0 

.  19 

.502 

-  .532 

5.6 

100.7 

24.0 

- 

.  14 

.459 

-  .344 

5.7 

103.0 

24.0 

- 

.  09 

.431 

-  .219 

5.8 

105.3 

24.0 

- 

.  6 

.414 

-  .137 

5.9 

107.6 

24.0 

- 

.  4 

.403 

-  .085 

6.0 

109.9 

24.0 

- 

.  3 

.396 

-  .  52 

V 
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TABLE 

B-2  (cont 

'd) 

(il)  Case  B, 

Mq=2.8,  Cooling 

V 

*1 

F1 

V 

°1 

Q  1 

W1 

0.0 

0.00000 

0.00000 

1.00882 

2.0000 

0.15119 

.2 

.  1002 

.20173 

1.00815 

2.0184 

.03311 

.4 

.  4000 

.40299 

1.00345 

2.0133 

-0.06425 

.6 

.  9030 

.60258 

0.99080 

1 .98493 

-  .19841 

.8 

.16182 

.79853 

.96649 

1.93441 

-  .30506 

1.0 

.25597 

.98818 

.92731 

1.86379 

-  .39348 

1.2 

.37458 

1.16831 

.87101 

1.77633 

-  .47226 

1.4 

.51980 

1.33539 

.79686 

1.67660 

-  .52037 

1.6 

.69391 

1.48594 

.70614 

1.57018 

-  .53857 

1.8 

.89894 

1.61697 

.60248 

1.46325 

-  .52570 

2.0 

1.13621 

1.72645 

.49181 

1.36180 

-  .48452 

2.2 

1.40583 

1 .81372 

.38163 

1.27089 

-  .42162 

2.4 

1.70625 

1.87965 

.27971 

1.19399 

-* .34620 

2.6 

2.0342 

1.92658 

.19253 

1 .13259 

-  .26810 

2.8 

2.385 

1.95788 

.12382 

1.08636 

-  .19577 

3.0 

2.7533 

1.97735 

.07412 

1.05352 

-  .13479 

3.2 

3.1342 

1.98862 

.  4119 

1.03153 

-  .08750 

3.4 

3.5232 

1.99467 

.  2121 

1.01764 

-  .  5357 

3.6 

3.9173 

1.99767 

.  1011 

1.00937 

-  .  3093 

3.8 

4.3143 

1.99905 

.  0447 

1.00472 

-  .  1684 

4.0 

4.7129 

1 .99964 

.  183 

1.00225 

-  .  0865 

4.2 

5.1122 

1.99987 

.  070 

1.00102 

-  .  420 

4.4 

5.5119 

1.99996 

.  25 

1.00044 

-  .  192 

4 .6 

5.9118 

1.99999 

08 

1.00018 

-  .  033 

4.8 

6.3118 

2.0000 

3 

1 .00007 

-  .  34 

5.0 

6.7117 

2.0000 

1 

1.00002 

-  .  13 

5.2 

7.1117 

2.0000 

0 

1 .00001 

-  .  05 

5.4 

7.5117 

2.0000 

0 

1.00000 

-  .  2 

5.6 

7.9H7 

2.0000 

0 

1.00000 

-  .  1 

5.8 

8.3117 

2.0000 

0 

1.00000 

-  .  0 

6.0 

8.7117 

2.0000 

0 

1.00000 

-  .  0 

O  OJ^i-VO  CO 
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TABLE  B-2  (cont'd) 


(ii)  Case  B,  Mn  =  2.8,  Cooling  (cont'd) 


V  f2  F2 

0.0000  0.0000 

.0843  1.7624 

.3540  3.825 

.8373  6.148 

1.5683  8.659 

1.0  2.589  11.255 

1.2  3.948  13.806 

1.4  5.700  16.166 

1.6  7.895  18.189 

1.8  10.568  19.756 

2.0  13-714  20.80 

2.2  17.275  21.31 

2.4  21.13  21.38 

2.6  ,25.11  21.13 

2.8  29.04  20.71 

5.0  32.80  20.26 

3.2  36.31  19.874 

3.4  39.56  19.586 

3.6  42.61  19.398 

3.8  45.49  19.288 

4.0  48.26  19.230 

4.2  50.97  19.203 

4.4  53.64  19.191 

4.6  56.29  19.186 

4.8  58.92  19.184 

5.0  61.56  19-184 

5.2  64.19  19.184 

5.4  66.82  19.184 

5.6  69.45  19.184 

5.8  72.08  19.184 

6.0  74.71  19.184 


V 

«2 

0  1 

W2 

8  .029 

0.0000 

0.7234 

9.585 

.1996 

1.1809 

11.010 

.4365 

1.1028 

12.156 

.6092 

0.5556 

12 .866 

.6361 

-  .3265 

12  .987 

.4697 

-1.3392 

12.404 

.1086 

-2.228 

11.073 

-  .3957 

-2.732 

9.057 

-  .9454 

-2.659 

6  .549 

-1.4163 

-1.9476 

3.862 

-1 .6893 

-0.7133 

1.3705 

-1.6846 

.7738 

-0.5824 

-1.3857 

2.172 

-1.7931 

-0.8413 

3.187 

-2.255 

-  .1466 

3.662 

-2.138 

.5888 

3.607 

-1.6984 

1.2712 

3.159 

-1.1745 

1.8402 

2.506 

-0  .7204 

2.273 

1.8198 

-  .3962 

2.575 

1.2182 

-  .1969 

2  .770 

0.7552 

-  .  .0890 

2  .887 

.4350 

-  .  367 

2.952 

.2334 

-  .  139 

2.986 

.1170 

-  .  049 

3.002 

.0549 

-  .  16 

3.010 

.  241 

-  .  05 

3.013 

.  100 

-  .  1 

3.014 

.  039 

-  .  0 

3-015 

.  14 

-  .  0 

3 .015 

.  05 

-  .  0 

3.015 

2 
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TABLE  B-2  (cont'd) 


(iii)  Case  C,  MQ=2.8,  Heating 


V 


1.0 

1.1 

1.2 
1.5 
1.4 


1.5 

1.6 

1.7 

1.8 

1.9 

2.0 

2.1 

2.2 

2.3 

2.4 

2.5 

2.6 

O  7 
*-  •  | 

2.8 

2.9 

3.0 

3.1 

3.2 

3.3 

3.4 


0.00000 
.  151 

.  608 
.  1379 
.  2473 


.  3903 
.  5680 
.  7820 
.10340 
.13259 


.16597 

.20378 

.24628 


.29373 


.34645 


.40473 

.46891 

.53934 

.61635 

.70030 


.79149 

.89023 

.99676 

1.11126 

1.23382 


1.36445 

1.50304 

1.64936 

1.80306 

1.96369 


2.1307 

2.3035 

2.4813 

2.6636 

2.8496 


F1 

V 

°1 

-0l* 

0.00000 

0.89885 

3.0000 

0.27834 

.  8988 

.89880 

2.9698 

.32520 

.17975 

.89849 

2.9350 

.37197 

.26956 

.89762 

2.8954 

.41848 

.35925 

.89592 

2.8513 

.46452 

.44871 

.89309 

2.8026 

.50981 

.53782 

.88885 

2.7494 

.55399 

.62642 

.88290 

2.6918 

.59667 

.71433 

.87494 

2  .6301 

.63737 

.80133 

.86471 

2.5644 

.67556 

.88719 

.85193 

2.4951 

.71068 

.97163 

.83635 

2.4224 

.74211 

1.05436 

.81778 

2.3468 

.76925 

1.13503 

.79602 

2.2687 

.79147 

1.21345 

.77098 

2.1887 

.80818 

1.28916 

.74260 

2.1073 

.81883 

1.36186 

.71089 

2.0251 

.82298 

1.43123 

.67598 

1.94292 

.82026 

1.49696 

.63806 

1.86133 

.81048 

1.55875 

.59745 

1.78106 

.79359 

1.61637 

.55454 

1.70284 

.76973 

1 .66960 

.50985 

1.62734 

.73925 

1.71830 

.46398 

1.55519 

.70271 

1.76238 

.41759 

1.48697 

.66085 

1.80183 

.37141 

1.42317 

.61459 

1.83669 

.32619 

1.36417 

.56500 

1.86712 

.28264 

1.31024 

.51324 

1 .89330 

.24144 

1.26155 

.46050 

1 .91550 

.20317 

I.21814 

.40798 

1.93405 

.16830 

1.17991 

.35677 

1.94929 

.13715 

1.14671 

.30788 

1.96161 

.10989 

1.11823 

.26211 

1.97139 

.08651 

1.09416 

.22010 

1.97903 

.  6689 

1.07408 

.18226 

1.98489 

.  5078 

1.05756 

.14881 
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TABLE  B-2  (cont'd) 


(iii)  Case  C,  Mq  =  2.8,  Heating  (cont'd) 


V 

fl 

F1 

V 

el 

-°1 

5.5 

5.0587 

1.98929 

0.03783 

1.04417 

0.11978 

5.6 

3.2504 

1.99255 

.  2765 

1.03546 

•09504 

5.7 

3.4242 

1.99490 

.  1983 

1.02503 

.  7433 

5.8 

3.6196 

1.99657 

.  1594 

1.01847 

.  5729 

5.9 

3.8162 

1.99774 

.  0961 

1.01346 

.  4352 

4.0 

4.0137 

1.99854 

.  650 

1.00968 

.  5258 

4.1 

4.2120 

1.99907 

.  431 

1.00686 

.  2403 

4.2 

4.4107 

1.99942 

.  280 

1.00480 

.  1747 

4.5 

4.6099 

1.99965 

.  178 

1.00332 

.  1251 

4.4 

4.8093 

1.99979 

.  Ill 

1 .00226 

.  0883 

4.5 

5.OO89 

1.99988 

.  068 

1.00152 

.  6i4 

4.6 

5.2087 

1.99995 

.  41 

1.00101 

.  421 

4.7 

5.4085 

1.99996 

24 

1.00066 

.  284 

4.8 

5 .6084 

1.99998 

14 

1.00042 

.  189 

4.9 

5.8083 

1.99999 

.  08 

1.00027 

.  124 

5.0 

6.0083 

2.0000 

4 

1.00017 

.  080 

5.1 

6.2083 

2.0000 

2 

1.00011 

.  51 

5.2 

6.4082 

2.0000 

1 

1 .00006 

.  52 

5.5 

6.6082 

2.0000 

1 

1 .00004 

20 

5.4 

6.8082 

2.0000 

0 

1 .00002 

12 

5.5 

7.0082 

2.0000 

0 

1.00001 

.  07 

5.6 

7.2082 

2.0000 

0 

1 .00001 

4 

5.7 

7.4082 

2.0000 

0 

1.00001 

2 

5.8 

7 .6082 

2.0000 

0 

1.00000 

1 

5.9 

7.8082 

2.0000 

0 

1.00000 

1 

6.0 

8.0082 

2.0000 

0 

1.00000 

0 

6.1 

8.2082 

2.0000 

0 

1.00000 

0 

6.2 

8.4082 

2.0000 

0 

1.00000 

0 

6.5 

8.6082 

2.0000 

0 

1.00000 

.  .  0 

6.4 

8.8082 

2.0000 

0 

1.00000 

0 

6.5 

9.0082 

2.0000 

0 

1.00000 

0 
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TABLE  B-2  (cont'd) 


(ill)  Case  C,  Mq=2.8,  Heating  (cont'd) 


V  f2  P 2 


0.0000 

0.0000 

.  131 

.7963 

.  547 

1.6747 

.1282 

2.634 

.2373 

3.670 

.3858  ■ 

4.781 

.5782 

5.961 

.8192 

7.202 

1.1137 

8.497 

1.4676 

9.836 

1.0  1.8870  11.207 

1.1  2.379  12.598 

1.2  2.951  13.993 

1.3  3.610  15.377 

1.4  4.367  16.735 

1.5  5.229  18.047 

1.6  6.207  19.298 

1.7  7.310  20.47 

1.8  8.547  21.55 

1.9  9.927  22.52 

2.0  11.455  23.36 

2.1  13.136  24.08 

2.2  14.972  24.66 

2.3  16.958  25.11 

2.4  19.086  25.42 

2.5  21.34  25.61 

2.6  23.71  25.69 

2.7  26.16  25.66 

2.8  28.68  25.56 

2.9  31.22  25.39 

3.0  33.76  25.19 

3.1  36.27  24.95 

3.2  38.73  24.71 

3.3  41.12  24.48 

3.4  43.43  24.26 


V 

q2 

V 

7.550 

0.0000 

0 .0679 

8.375 

.  198 

.3057 

9.190 

.  569 

.4161 

9.984 

.  988 

.4025 

10.746 

.1335 

.2706 

11.462 

.1493 

.0285 

12.117 

.1358 

-  .3128 

12.697 

.0839 

-  .7387 

13.187 

-  .  141 

-1.2316 

13-570 

-  .1639 

-1.7699 

13.830 

-  .3688 

-2.329 

13.954 

-  .6295 

-2.881 

13.927 

-  .9439 

-3.397 

13.738 

-1.3067 

-3.845 

13.379 

-1.7096 

-4.194 

12.847 

-2.141 

-4.415 

12.142 

-2.587 

-4.482 

11.271 

-3.032 

-4.375 

10.248 

-3.456 

-4.081 

9.093 

-3.842 

-3.598 

7.834 

-4.170 

-2.932 

6.504 

-4.423 

-2.102 

5.142 

-4.586 

-1.1383 

3.791 

-4.647 

-0.0793 

2.495 

-4.600 

1.0275 

1.2970 

-4.442 

2.130 

0.2333 

-4.176 

3.176 

-  .6654 

-3.810 

4.117 

-1.3796 

-3.357 

4.910 

-1.9011 

-2 .834 

5.523 

-2.234 

-2.259 

5.938 

-2.393 

-1.6531 

6.149 

-2.403 

-1.0359 

6.163 

-2.294 

-0.4263 

6.001 

-2.099 

.1594 

5.690 
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TABLE  B-2  (cont'd) 

( Hi)  Case  C,  Mq  =  2.8,  Heating  (cont'd) 


V 

f2 

F2 

V 

0 

2 

3.5 

45.65 

24.06 

-1.8495 

0.7078 

3 .6 

47.78 

23.89 

-1.5756 

1.2092 

3.7 

49.85 

23.74 

-1.3005 

1.6572 

3.8 

51.79 

23.63 

-1.0420 

2  .049 

3-9 

55.68 

23.53 

-0  .8116 

2.384 

4.0 

55.51 

23.46 

-  .6151 

2.666 

4.1 

57.29 

23.41 

-  .4540 

2.897 

4.2 

59.02 

23.37 

-  .3266 

3.084 

4 .3 

60.72 

23.34 

-  .2291 

3.232 

4.4 

62.59 

23.32 

-  .1568 

3.347 

4.5 

64.04 

23.31 

-  .1047 

3.434 

4.6 

65  .68 

23.30 

-  .0683 

3.500 

4.7 

67.50 

23.30 

-  •  435 

3.549 

4 .8 

68.92 

23.29 

-  .  270 

3.584 

4.9 

70.55 

23.29 

-  .  164 

3.609 

5.0 

72.15 

23.29 

-  •  097 

3 .626 

5 .1 

75.74 

23.29 

-  .  56 

3.638 

5.2 

75.54 

23.29 

-  •  31 

3.646 

5-3 

76.95 

23.29 

-  •  17 

3.652 

5.4 

78.55 

23.29 

-  .  08 

3.655 

5.5 

80.15 

23.29 

-  .  4 

3 .658 

5.6 

81.75 

23.29 

-  .  1 

3.659 

5.7 

85.52 

23.29 

0 

3.660 

5.8 

84.92 

23.29 

1 

3 .661 

5.9 

86.52 

23.29 

1 

3 .661 

6.0 

88.11 

23.29 

1 

3 .661 

6.1 

89.71 

23.29 

1 

3  .661 

6.2 

91.51 

23.29 

A 

•  a 

3 .661 

6.3 

92.90 

23.29 

2 

3  .661 

6.4 

94.50 

23.29 

2 

3 .662 

6.5 

96.10 

23.29 

2 

3 .662 

5  .262 
4.754 
4.201 
3.653 
3.078 

2.556 
2.081 
1.6650 
1.5041 
1 .0041 

0.7595 

.5641 

.4118 

.2955 

.2084 

.1446 

.0987 
.  662 
.  458 

.  285 

.  185 
.  116 
.  075 
.  45 

.  28 

.  18 
.  11 
.  08 
.  6 
.  4 

.  4 
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TABLE  B-2  (cont‘d) 


V 


1.0 

1.1 

1.2 

1.3 

1.4 


1.5 

1.6 

1.7 

1.8 
1.9 


2.0 

2.1 

2.2 

2.3 

2.4 


2.5 

2.6 

2.7 

2.8 

2.9 

3.0 

3-1 

3.2 

3.3 

3.4 


(Hi) 

Case  C, 

MQ  =  2.8, 

Heating  (cont’d) 

f3 

p? 

V 

®3 

e3* 

0.000 

-  .053 

-  .182 

-  .345 

-  .496 

0.00 

-2.90 

-4.41 

-4.51 

-3.14 

-35.8 
-22.1 
-  8.10 
6.30 
21.0 

0.000 

.618 

1.553 

2.77 

4.21 

4.54 

7.80 

10.83 

13.39 

15.23 

-  .586 

-  .565 

-  .377 
.040 
.752 

-0.293 

4.07 

9.95 

17-35 

26.3 

36.0 

51.2 

66.4 

81.6 

96.6 

5.79 

7.41 

8.97 

10.36 

11.48 

16.18 

16.11 

14.94 

12.69 

9.50 

1.855 

3.37 

5.45 

8.18 

11.66 

36 .6 

48.5 

61.6 
76.0 
91.4 

111.1 

125.0 

137.8 

149.4 

159.0 

12.24 

12.58 

12.48 

11.97 

11.12 

5.56 

1.218 

-3.12 

-6.95 

-9.76 

16.01 

21.3 

27.7 

35-3 

44.1 

107.7 

124.6 

141.6 

158.4 

174.4 

166.2 

170.2 

170.0 

164.8 

153.4 

10.07 

8.99 

8.10 

7.65 

7.88 

-11.02 

-10.25 

-  7.12 

-  1.472 

6 .63 

54.1 

65.2 

77.3 

90.1 

103.1 

188.8 

201 

210 

216 

217 

135.0 

108.9 

75.1 

34.3 

-11.47 

9.04 

11.30 

14.77 

19.45 

25.3 

16.84 

28.5 

40.8 
52.7 

62.9 

115.9 

127.9 
138.5 
147.2 

153 .4 

214 

205 
192.5 
176.3 
157  .5 

-  59.4 
-105.8 
-146.3 
-177.4 
-196.2 

32.0 

39.2 

46.7 

53.9 

60.3 

70.5 
74.3 
73-9 
68.9 

59.6 

157.0 

157.9 

156.3 

152.6 

147.1 

137.5 

117.6 
99.0 

82.5 

68.6 

-  20.2 
-194.3 
-176.7 
-152.2 
-124.4 

65.7 

69-5 

71.8 

72.2 

70.9 

46.6 
30.8 
13.38 
-  4.28 
-21.0 
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TABLE  B-2  (cont'd) 


(iii)  Case  Cj  Mq-2.8 ,  Heating  (cont'd) 


9  f3 

f3 

V 

e3 

6  ' 

3 

3.5  140.5 

3.6  153.5 

3.7  126.4 

3.8  119.7 

3.9  113.8 

57.6 

49.2 
43.1 
39.0 

36.3 

-96.7 

-71.5 

-50.3 

-33.6 

-21.3 

68.1 

63.9 

58.7 

52.8 

46  .5 

-35.6 

-47.5 

-56.2 

-61.4 

-63.5 

4.0  108.8 

4\1  104.9 

4.2  102.1 

4.3  100.2 

4.4  99.3 

34.6 

33.6 
33.0 
32.8 

32.6 

-12.84 

-  7.31 

-  3.96 

-  2.06 

-  1.088 

40.2 

34.1 

28.3 

23.2 
18.64 

-62.6 

-59.5 

-54.6 

-48.7 

-42.2 

4.5  99.3 

4.6  99.9 

4.7  101.1 

4.8  102.8 

4.9  104.9 

32.5 

32.5 

32.4 

32.4 

32.4 

-  0.636 

-  .451 

-  .384 

-  .354 

-  .326 

14.75 

11.51 

8.87 

6.77 

5.13 

-35.6 

-29.3 

-23.6 

-18.56 

-14.29 

5.0  107.2 

5.1  109.8 

5.2  112.5 

5.3  115.3 

5.4  118.2 

32.3 

32.3 

32.3 

32.3 

32.2 

-  .290 

-  .246 

-  .199 

-  .154 

-  .115 

3.89 

2.95 

2.27 

1.783 

1.439 

-10.78 

-  7.96 

-  5.77 

-  4.10 

-  2.86 

5.5  121.2 

5.6  124.2 

5.7  127.2 

5.8  130.2 

5.9  133.3 

32.2 

32.2 

32.2 

32.2 

32.2 

-  .083 

-  .  57 

-  .  39 

-  .  26 
-  .  16 

1.201 

1.040 

0.932 

.862 

.816 

-  1.954 

-  1.312 

-  0 . 866 

-  .561 

-  .358 

6.0  136.3 

6.1  139.4 

6.2  142.5 

6.3  145.6 

6.4  148 .6 

32.2 

32.2 

32.2 

32.2 

32.2 

-  .  10 
-  .  06 

-  .  4 

-  .  2 

-  .  2 

.788 

•  770 

•  759 
.753 
.749 

-  .224 

-  .138 

-  .084 

-  .  51 

-  .  31 

6.5  151.7 

32.2 

-  .  1 

.746 

-  .  19 

00O\4=T0  o 
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TABLE  B-2  (cont'd) 


(iv)  Case  D,  Mq=2,  Adiabatic 

?  f  1  Fx  F  '  ©1  -©! ' 

0.00000  0.00000  1.09758  1.68094  0.00000 

.  1309  .21947  1.09663  1.67381  .07125 

.  5265  .43827  1.08993  1.65250  .14159 

.11954  .65468  1.07182  1.61739  .20867 

.21505  .86588  1.03706  1.56949  .26878 


1.0 

.34079 

1 .06811 

0.98150 

1.51065 

.31726  ; 

1.2 

.49843 

1.25694 

.90299 ■ 

1.44368 

.34934 

1.4 

.68947 

1.42781 

.80230 

1.37226 

.36133 

1.6 

.91476 

1.57666 

.68378 

1.30059 

.35183  j 

1.8 

1.17407 

1 .70064 

•55520 

1.23285 

.32255  1 

2.0 

1.46570 

1 .79873 

.42661 

1.17257 

.27827  1 

2.2 

1.78632 

1.87196 

.30833 

1.12209 

.22583  ! 

2.4 

2.1313 

1.92328 

.20849 

1.08231 

.17239 

2.6 

2.4953 

1.95686 

.13134 

1.05280 

.12384 

2.8 

2.8732 

1.97730 

.07684 

1.03221 

.08376 

3.0 

3.2603 

1.98885 

.  4166 

1.01886 

.  5337 

3.2 

3.6532 

1.99489 

.  2092 

1.01026 

.  3206 

3-4 

4.0495 

1.99782 

.  0972 

1.00535 

.  1815 

3-6 

4.4477 

1.99913 

.  419 

1.00265 

.  0969 

3.8 

4.8468 

1 .99968 

.  167 

1.00124 

.  488 

4.0 

5.2464 

1.99989 

.  062 

1.00055 

.  232 

4.2 

5.6463 

1.99996 

22 

1.00023 

.  104 

4.4 

6.0462 

1.99999 

.  07 

1.00009 

.  044 

4.6 

6.4462 

2.0000 

2 

1.00004 

18 

4.8 

6.8462 

2.0000 

1 

1.00001 

.  07 

5.0 

7.2462 

2.0000 

0 

1.00000 

2 

5-2 

7.6462 

2.0000 

0 

1.00000 

1 

5.4 

8.0462 

2  .0000 

0 

1 .00000 

0  r: 

5.6 

8.4462 

2.0000 

0 

1.00000 

0 

5.8 

8.8462 

2.0000 

0 

1.00000 

o 

i  ;• 

6.0 

9.2462 

2  .0000 

0 

1.00000 

0  !! 

1 
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TABLE  B-2  (cont’d) 


(iv)  Case  D,  M0=2,  Adiabatic  (cont'd) 


V 

f2 

P2 

V 

02 

Q  » 
w2 

0.0 

0.0000 

0.0000 

4.034 

0.6071 

0.0000 

.2 

.  524 

1.0403 

6.356 

.6865 

.7200 

.4 

.2443 

2.532 

8.519 

.8661 

1.0065 

.6 

.6264 

4.425 

10.330 

1.0630 

0.9046 

.8 

1.2490 

6 .626 

11.572 

1.2083 

.5104 

1.0 

2.162 

9.002 

12.047 

1.2576 

-  .0253 

1.2 

3.412 

11 .384 

11 .617 

1.2011 

-  .5165 

1.4 

5.036 

13.587 

10.270 

1.0660 

-  .7854 

1.6 

7.049 

15.441 

8.162 

0.9095 

-  .7210 

1.8 

9.432 

16.823 

5.618 

.8001 

-  .3241 

2.0 

12.121 

17.687 

3.070 

.7944 

.2886 

2.2 

15.011 

18.077 

0.9319 

.9176 

.9323 

2.4 

17.978 

18.104 

-  .5257 

1.1576 

1.4309 

2.6 

20.91 

17.914 

-1.2514 

1.4735 

1.6815 

2.8 

23.73 

17.641 

-1.3879 

1.8134 

1.6745 

3-0 

26.40 

17.379 

-1.1732 

2.131 

1.4722 

3.2 

28.92 

17.177 

-0.8329 

2.396 

1.1669 

3.4 

31.34 

17.043 

-  .5154 

2.597 

0.8441 

3.6 

33.66 

16.964 

-  .2832 

2.737 

.5616 

3.8 

35.93 

16.923 

-  .1398 

2.827 

.3454 

4.0 

38.17 

16.903 

-  .0625 

2  .880 

.1971 

4.2 

40.39 

16.895 

-  .  255 

2.910 

.1047 

4.4 

42.61 

16.892 

-  .  095 

2.925 

.0519 

4.6 

44.81 

16.891 

-  .  33 

2.932 

.  241 

4.8 

47.02 

16.890 

-  .  11 

2.935 

.  104 

5.0 

49.22 

16.890 

-  .  03 

2.937 

.  043 

5.2 

51.42 

16.890 

-  .  1 

2.937 

.  16 

5.4 

53.63 

16.890 

-  .  0 

2.937 

.  06 

5.6 

55.83 

16.890 

-  .  0 

2.937 

2 

5.8 

58.03 

16.890 

-  .  0 

2.937 

1 

6.0 

60.24 

16.890 

0 

2.937 

0 

37 


TABLE  B-2  -(cont'd) 


(v)  Case  E,  Mq=2,  Cooling 


V 

h 

pi 

V 

el 

0  ' 
W1 

0.0 

.2 

.4 

.6 

.8 

0.00000 
.  1568 
.  6248 
.14058 
.25074 

0.00000 

.22853 

.45624 

.68099 

.89941 

1.14297 

1.14177 

1.13347 

1.11127 

1.06921 

1.44500 

1.46182 

1.46320 

1.44949 

1.42167 

0.12277 

.  4542 

-  .  3131 

-  .10500 

-  .17156 

1.0 

1.2 

1.4 

1.6 

1.8 

.39403 

.57157 

.78408 

1.03144 

1.31224 

1.10706 

1.29890 

1.46998 

1 .61623 

1.73526 

1.00301 

0.91122 

.79616 

.66437 

.52585 

1.38168 

1.33249 

1.27796 

1.22238 

1.16984 

-  .22584 

-  .26280 

-  .27890 

-  .27345 

-  .24915 

2.0 

2.2 

2.4 

2.6 

2.8 

1.62354 

1.96097 

2.3193 

2.6931 

3.0776 

1 .82689 

1 .89320 
1.93808 
1.96637 

1 .98294 

.39227 

.27423 

.17882 

.10839 

.06093 

1.12360 

1.08559 

1.05634 

1.03523 

1.02092 

-  .21166 

-  .16805 

-  .12499 

-  .08725 

-  •  5726 

3.0 

3.2 

3.4 

3.6 

3.8 

3.4690 

3.8644 

4.2621 

4.6610 

5.0605 

1.99192 

1 .99644 

1 .99854 
1.99944 
1.99980 

.  3172 
.  1528 
.  0682 
.  282 
.  108 

1.01179 

1.00631 

1 .00320 
1.00154 
1.00070 

-  .  3538 

-  .  2060 

-  .  1131 

-  .  0585 

-  .  286 

4.0 

4.2 

4.4 

4.6 

4.8 

5.4603 

5.8602 

6.2602 

6.6602 

7.0602 

1.99993 

1.99998 

1 .99999 
2.0000 

2 .0000 

.  039 

.  13 

04 

1 

0 

1.00030 

1.00012 

1.00005 

1.00002 

1.00001 

-  .  132 

-  .  057 

-  .  24 

-  .  09 

-  .  3 

3.0 

5.2 

5.4 

5.6 

5.8 

7.4602 

7.8602 

8.2602 

8.6602 

9.0602 

2 .0000 

2 .0000 
2.0000 
2.0000 
2.0000 

0 

0 

0 

0 

0 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

-  .  1 

-  .  0 

-  .  0 

-  .  0 

-  .  0 

6.0 

9.4602 

2.0000 

0 

1.00000 

-  .  0 

Oocr»-»=-ro  o 
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TABLE  B-2  (cont'd) 


V 


1.0 

1.2 

1.4 
1.6 
1.8 

2.0 

2.2 

2.4 
2.6 
2.8 

3.0 

3.2 

3.4 

3.6 

3.8 

4.0 

4.2 

4.4 

4.6 

4.8 

5.0 

5.2 

5.4 

5.6 

5.8 

6.0 


(V) 

Case  E, 

M0=2,  Cooling 

(cont'd) 

f2 

F2 

V 

02 

V 

0.0000 

.0712 

.3173 

.7838 

1.5138 

0.0000 

1.1284 

2.689 

4.622 

6.823 

4.514 

6.754 

8.805 

10.442 

11.437 

0.0000 

.1690 

.4575 

.7199 

.9462 

0.4449 

1.1688 

1.4441 

1.3206 

0.9055 

2.548 

3.921 

5.651 

7.732 

10.117 

9.141 

11.396 

15-401 

15-000 

16.099 

11.590 

10.800 

9.120 

6.789 

4.207 

1.0735 

1.0956 

1.0431 

0.9723 

.9445 

.3628 

-  .1152 

-  .3613 

-  .2934 

.0525 

12.724 

15.442 

18.161 

20.80 

23.30 

16.695 
16 .867 
16.750 
16.490 
16.204 

1 .8269 

0.0176 

-1.0519 

-1.4570 

-1.3473 

1.0029  > 
1.1588 

1 .3925 
1.6653 
1.9362 

.5407 

1.0009 

1.3011 

1.3895 

1.2902 

25.68 

27.93 

30.10 

32.21 

54.28 

15.964 

15.792 

15.686 

15.628 

15.599 

-1.0299 

-0.6775 

-  .3930 

-  .2039 

-  .0955 

2.174 

2.362 

2.499 

2.589 

2.645 

1.07H 

0.8076 

.5585 

.3563 

.2106 

36.34 

38.38 

30.42 

42.45 

44.48 

15.586 

15.581 

15.579 

15.578 

15.578 

-  .  407 

-  .  158 

-  .  057 

-  .  19 

-  .  06 

2.677 

2.694 

2.702 

2.706 

2.708 

.1158 
.0593 
.  284 
.  127 
.  053 

46.51 

48.55 

50.58 

52.61 

54.64 

15.578 

15.578 

15.578 

15.578 

15.578 

-  .  2 

-  .  0 

-  .  0 

-  .  0 

-  .  0 

2.709 

2.709 

2.709 

2.709 

2.709 

.  21 
.  08 

.  3 

1 
0 

56 .67 

15.578 

-  .  0 

2.709 

0 
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TABLE  B-2  (cont'd) 


(vi)  Case  F,  MQ=2,  Heating 


V 


0.0 

0 .00000 

0.00000 

.2 

.  957 

.20422 

.4 

.  5900 

.40798 

.6 

.  8962 

.61008 

.8 

.16311 

.80853 

1.0 

.26147 

1.00055 

1.2 

.38691 

1.18269 

1.4 

.54174 

1.55123 

1.6 

.72803 

1.50243 

1.8 

.94729 

1.63515 

2.0 

1.19994 

1.74130 

2.2 

1.48495 

1.82643 

2.4 

1.79957 

1.88975 

2.6 

2.1396 

1 .93400 

2.8 

2.4998 

1.96290 

5.0 

2.8749 

1.98048 

4.2 

3.2602 

1.99041 

5.4 

5.6521 

1.99561 

5.6 

4.0478 

1.99813 

5.8 

4.4457 

1.99925 

4.0 

4.8447 

1.99972 

4.2 

5.2442 

1.99990 

4.4 

5.6440 

1.99997 

4.6 

6.0439 

1.99999 

4.8 

6.4439 

2.0000 

5.0 

6.8439 

2.0000 

5.2 

7.2439 

2.0000 

5.4 

7.6439 

2.0000 

5.6 

8.0439 

2.0000 

5.8 

8.4439 

2.0000 

6.0 

8.8439 

2.0000 

V 

el 

-Q  * 
W1 

1.02125 

2,1670 

0.22692 

1 .02060 

2.1154 

. 28852 

1.01600 

2.0517 

.34888 

1.00338 

1.97612 

.40580 

0.97873 

1.88979 

.45610 

.95845 

1.79439 

.49578 

.87995 

1.69247 

.52060 

.80232 

1.58740 

.52675 

.70703 

1.48317 

.51191 

.59831 

1.38404 

.47610 

.48297 

1.29394 

.42226 

.36945 

1.21597 

.35598 

.26615 

1.15190 

.28449 

.17958 

1.10203 

.21511 

.11299 

1.06533 

.15367 

.06606 

1.03980 

.10364 

.  3582 

1.02305 

.06596 

•  1799 

1.01268 

.  3960 

.  0836 

1.00661 

.  2242 

.  360 

1.00327 

.  1198 

.  144 

1.00153 

.  0604 

.  054 

1.00068 

.  287 

.  19 

1.00029 

.  129 

06 

1.00011 

•  055 

2 

1.00004 

22 

1 

1 .00002 

08 

0 

1.00001 

3 

0 

1.00000 

1 

0 

1.00000 

0 

0 

1.00000 

0 

0 

1.00000 

0 
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TABLE  B-2  (cont'd) 


(vi)  Case  F,  Mq  =  2,  Heating  (cont'd) 


V 


0.0 

0.0000 

0.0000 

.2 

.  378 

.8969 

.4 

.1830 

2.280 

.6 

.4845 

4.109 

.8 

.9945 

6.311 

1.0 

1 .7687 

8.772 

1.2 

2.866 

11.346 

1.4 

4.343 

13.857 

1.6 

6.249 

16.127 

1.8 

8.605 

17.998 

2.0 

11.396 

19.364 

2.2 

14.551 

20.20 

2.4 

17.953 

20.56 

2.6 

21.45 

20.57 

2.8 

24.91 

20.37 

3-0 

28.23 

20.09 

3.2 

31.36 

19.835 

3.4 

34.31 

19.636 

3.6 

37.11 

19.505 

3.8 

39.80 

19.429 

4.0 

42.42 

19.389 

4.2 

44.99 

19.370 

4.4 

47.54 

19.362 

4.6 

50.08 

19.359 

4.8 

52.61 

19.358 

5.0 

55.13 

19.358 

5.2 

57.66 

19.358 

5-4 

60.18 

19.358 

5.6 

62.71 

19.358 

5.8 

65.23 

19.358 

6.0 

67.76 

19.358 

V 

02 

0  1 

2 

3.239 

0.0000 

0.3602 

5.721 

.1452 

1.0168 

8.078 

.3775 

1.2360 

10.154 

.6133 

1 .0621 

11.766 

.7818 

0.5829 

12.724 

.8351 

-  .0636 

12  .865 

.7571 

-  .6978 

12.100 

.5697 

-1.1278 

10.465 

.3303 

-l.i979 

8.151 

.1193 

-0.8423 

5.496 

.0182 

-  .1200 

2.919 

.  843 

.7943 

0 .8084 

.3333 

1.6686 

-  .5999 

.7355 

2.297 

-1.2796 

1.2288 

2.570 

-1.3855 

1.7405 

2.490 

-1.1570 

2.209 

2.153 

-0.8147 

2.594 

1.6873 

-  .5010 

2.884 

1.2102 

-  .2738 

3.084 

0.7996 

-  .1345 

3.211 

.4889 

-  .0598 

3.286 

.2776 

-  .  243 

3.328 

.1468 

-  .  090 

3.349 

.0725 

-  •  31 

3-359 

.  335 

-  .  10 

3.363 

.  145 

-  .  03 

3.305 

.  059 

-  .  1 

3.366 

.  23 

-  .  0 

3.366 

.  08 

-  .  0 

3.366 

.  3 

-  .  0 

3.367 

1 
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TABLE  B-3 


Values  o ? /3n>  tn,  and  kR 


Case 

A 

/32 

A 

tl 

t2 

A 

11.3357 

47.20 

198.4 

3.2388 

0.373 

32.4 

B 

10.5533 

38.29 

3.0152 

-0.424 

C 

12.8111 

61.16 

310 

3  .6603 

O.727 

51.2 

D 

10.2810 

57.72 

2.9374 

5.71 

E 

9.4821 

46.50 

2.7092 

4.11 

P 

11.7828 

76.93 

3.366 5 

7.81 

Case 

*1 

<*2 

kl 

k2 

k3 

A 

8.0969 

20.60 

96.2 

-1.03277 

-0.652 

-11.02 

B 

7.5380 

15.986 

-0.96149 

-0.327 

C 

9.1508 

26.94 

153.4 

-1.16719 

-0.913 

-17.41 

D 

7.3435 

30.45 

-1.83589 
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FIG.  2 -VARIATION  OF  BOUNDARY-LAYER  THICKNESS, 
DISPLACEMENT  THICKNESS,  AND  MOMENTUM 
THICKNESS  ALONG  TUBE  LENGTH 


HEAT-TRANSFER  RATE,  qD/X0T0  SHEAR-STRESS  RATIO,  r/r 


RELATIVE  DISTANCE  RATIO,  r, 


FIG. 4- PROFILES  OF  STAGNATION  TEMPERATURE, 
SHEAR  STRESS,  AND  HEAT -TRANSFER  RATE 
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FIG.  6  -VARIATION  OF  RECOVERY  FACTOR,  FRICTION 
COEFFICIENT,  MEAN  PRESSURE  -  GRADI  ENT 
COEFFICIENT,  AND  HEAT-TRANSFER  COEFFICIENT 
ALONG  TUBE  LENGTH 
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